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ABSTRACT 

The formation channel of the tens of compact debris discs which orbit white dwarfs (WDs) 
at a distance of one Solar radius remains unknown. Asteroids that survive the giant branch 
stellar phases beyond a few au are assumed to be dynamically thrust towards the WD and 
tidally disrupted within its Roche radius, generating extremely eccentric (e > 0.98) rings. 
Here, we establish that WD radiation compresses and circularizes the orbits of super-micron 
to cm-sized ring constituents to entirely within the WD’s Roche radius. We derive a closed 
algebraic formula which wcll-approximates the shrinking time as a function of WD cooling 
age, the physical properties of the star and the physical and orbital properties of the ring 
particles. The shrinking timescale increases with both particle size and cooling age, yielding 
age-dependent WD debris disc size distributions. 

Key words: minor planets, asteroids: general stars: white dwarfs - methods: numerical - 
celestial mechanics - planet and satellites: dynamical evolution and stability protoplanetary 
discs 


1 INTRODUCTION 


Mounting discoveries of debris orbiting white dwarfs (WDs) 
presage increased scrutiny of post-ma in-sequence planetary 
The over 30 
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3arber et all 2014: Bergfors et all 

2014; 
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200G 

200?l. 20081: Gansicke 201 ll: Farihi et all 20121; Melis et al. 

2012; 

Wilson et all 20141: Manser et al. 2015f) so far detected all have 


radial extents of just about 1 Solar radius ~ 0.005 au < 10° 
km. Such compact configurations are absent from main sequence 
planetary studies because they cannot exist; the discs would be 
inside of the star! Consequently, achieving an understanding of 
how these post-main-sequence discs are formed requires a dif¬ 
ferent approach. 

Protoplanetary discs orbiting young main-sequence stars 
form out of a collapsing stellar birth cloud, but WD discs 
must instead be formed fro m the tidal disruption of objects 
which encounter t he WD ([Gra ham et all Il990l : ljur 2 120031: 
iDebes et al.l 120121 : iBear fc Sokerl 2013fh However, because of 
mass loss and stellar tides on the giant branch phases of stel¬ 
lar evolution, there exists a “planet desert”. Mass loss serves 
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to push indi vidual planets ou t ward (sometimes to the interstel- 
lar m edium; yeraseUahl kOlll : IVeras fc Tout|[2012l : lAdams et al.l 
120131 : IVeras et a l. 201da) , and may cause multiple planets to be¬ 
come unstable (|DunciniJ^J(<issangj Il998|; iDebe^&Sigurdsson 


2002; Veras_et_alJ_2013aP Vovatzis et al] 20131 : Mustill et al 


2014 ; IVeras fc Gansickell2015l '). leading to collisions, engulfment 


within the star, or escape from the system. Concurrently, tidal 
effects between the giant star and a close-in planet could over¬ 
come the outward evolu tion due to mass loss and cause the star 
to swallow the planet (fVillaver fc Livio 2009; Kunitomoet_aH 


2011 1 _ Mustill_&Villayer 2012; Spiegel fc Madhusudhan 20121 ; 


Adams fc Bloc.hl 120131: IVillaver et al.ll20ljh~ Therefore, we 

pect a void of planetary bodies within a few au of WDs. 

Asteroids, comets, moons and planets are instead likely to 
exist beyond a few au in WD systems. Remnant exo-asteroid or 
exo-Kuiper belts can then dynamically interact with these plan¬ 
ets, flinging astero i ds towards the WD and eventual destruc- 
tion ( Bonsor et al] 1201 ll ; IDebes et al.l 120121 ; iFrewen fc Hansenl 
1201# . The d sruption of a single asteroid, which must be on 
a highly eccentric orbit in order to enter the disruption sphere 


1 Other potential reservoirs of d isrupted material, such as exo-Oort 
clouds , are both compositionally ([Z uckerman et al. 2007|) and dynam¬ 
ically llVeras et alj|2014bl ; I Stone et al .Il2015h unlikely to represent the 
primary source of the discs. 
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or Roche radius) o f the WD, results in a highly eccentric ring 
Veras et al. 120143 . hereafter Paper I). By modelling the aster¬ 
oid as a rubble pile composed of equal-mass, equal-radii inde¬ 
structible spheres, Paper I found that the resulting eccentric ring 
is collisionless. The tidal disruption of multiple asteroids in quick 
succession might instead produce collisional rings or discs. How¬ 
ever, here we treat the single asteroid (collisionless) case as the 
next logical step to follow Paper I. The resulting eccentric ring 
maintains the same orbit as the progenitor, but precesses due to 
general relativity. Paper I found that the eccentricity e of these 
orbits always satisfies e > 0.98, with typical values of e ~ 0.995; 
an asteroid with a semimajor axis a of 10 au that skims a typical 
WD would satisfy e ~ 0.99999. If the rings remain collisionless, 
then without the influence of additional forces, they will never 
change shape. 

However, observations of gaseous WD discs reveal circular 
or near-circular geometries. Double-peaked metal line emission 
profiles suggest e ~ 0 .02 and e ~ 0.2 in th e two most robustly- 
constrained systems dGansicke et al.l 120061 . l2008l l. The process 
of converting a highly eccentric orbit into a nearly circular one, 
while drastically reducing the semimajor axis, has rarely been 
discussed in the context of WDs, and is the issue we address in 
this paper. 

To do so, we introduce radiative forces into the two-body 
problem, as each ring constituent can be considered to be orbit¬ 
ing the WD unperturbed by any other body if the constituents 
are indeed collisionless. We focus on radiation and do not con¬ 
sider effects from tidal dissipation in this work. Regardless, radi¬ 
ation forces may dominate over tidal dissipation effects because 
the latter are likely to cause orbital changes of a ring particle 
on a much longer timescale. The effects of tidal dissipation in 
small and light ring constituents should be negligible, as can 
be d educed directly from the e quations of motion (e.g. equation 
8 of lBeauee fo Nesvornvl 1201^ 1 or from relations for the orbit- 
averaged semimajor axi s and eccentricity dam ping timescales 
(e.g. equations 10-11 of iNagasawa fc Ida 120111) . Both of those 
studies account for dynamical tides, a physical description in 
which highly eccentric orbits may be treated. 

Radiation, in the restricted, specific case of Poynting- 
Robertson drag, has become a crucial component of explain- 
ing accretion from compac t circum s tellar dis cs onto WDs 
llBochkarev fe Rafikovl 1201 ll : iRafikovI l2011al lbl: iMetzger et al.l 
l2012h . [Wyatt et al. FErnl)' has constrained the possible size dis¬ 
tributions of the accreting material by fitting model param¬ 
eters to observationally-inferred accretion rate distributions. 
We know that accretion onto the WD occurs because associ¬ 
ated with every dusty and gaseous disc are unmistakable sig¬ 
natures of photospheric metal pollution. These metals cannot 
arise from stellar dredge-up nor from the interstellar medium 
(ISM) due to the high accretion rate, the distribution of inter¬ 
stellar c louds and the existence of metal-rich WDs that lack hy¬ 
drogen dA annest ad et al. 1993|_ Friedrich_etaL 20041: Jural 200d ; 


iKilic fc Redfieldll2007l : iGansicke et alj|2008l : iFarihi et alTl2010T > . 

Here, we consider WD radiation in a different, earlier context, 
before a circular compact disc has formed. We are particu¬ 
larly motivated by the suggestion that Poynting-Robertson drag 
plays an important role in the dynamics of highly-eccentric or¬ 
bits (IVokrouhlickv et all 1200 il l. as the tidally-disrupted debris 
from Paper I settled into such orbits. 

In Section 2, we present our framework of perturbations 


due to WD radiation. Sections 3 and 4 describe how radiation 
changes orbits on short and long timescales, respectively. Based 
on these results, we derive the key formulae for orbital shrinkage 
in Section 5 before discussing the implications for the resulting 
ring structure due to collisions and the Yarkovsky effect in Sec¬ 
tion 6 and concluding in Section 7. 


2 GENERAL FORMULATION 


We let f(t) represent the distance from the centre of the ring 
particle to the centre of the WD, v(t) the velocity of the par¬ 
ticle with respect to the centre of the WD, m the particle’s 
mass, A the particle’s momentum-carrying cross-sectional area, 
c the speed of light and L(t) the time-varying luminosity of the 
WD. Then t he extra accelerati on of the particle due to stellar 
radiation is (IVeras et al.ll2015al 'l 



AL 

4nmcr 2 


Qabsl + Q ref H + fcQyarY 




(1) 


where Q a bs is the absorption efficiency, Q re f is the reflection ef¬ 
ficiency, Q yar = Qabs — Qref, I is the 3x3 unit matr ix, and ¥ is 
the 3x 3 Yarkovsky matrix (given in equation 28 of IVeras et al.l 
l2015al b The thermal redistribution parameter k satisfies 0 ^ 
k ^ 1/4. The three terms in equation Q refer to the contri¬ 
butions from radiation which is absorbed, directly reflected and 
thermally emitted. The vector 



( 2 ) 


represents the relativistically-corrected direction of incoming ra¬ 
diation. 

The above equations hold when the particle is at least a 
few orders of magnitude larger than the wavelength of incoming 
radiation, A. For WDs, A changes with cooling age (time since 
the star became a WD). Young (~4 Myr old), very hot WDs 
with temperatures of about 40000 K will have a value of A that 
is just under 1 x 10 — ' m. Alternatively, old (~2 Gyr old) and 
cold (6000 K) WDs have A ~ 5x 10 -7 m. Effectively, throughout 
a WD’s lifetime, A varies by less than an order of magnitude, 
and is always safely smaller than a micron. 

The contribution from the Yarkovsky term is negligible 
when the particle size is larger than about 1 cm - 10 m 
(IVeras et al.ll2015ah and when the particle spins, such that some 
of the absorbed radiation is redistributed before being emit¬ 
ted. When the Yarkovsky effect is active, it can represent the 
dominant contribution to the motion: the resulting perturbation 
may eject an asteroid orbiting a 10 3 Lq evolved star (such as an 
asymptotic giant branch star or very young WD) in just tens of 
million years. Here, we assume that the tidal disruption of an 
asteroid into sub-metre-sized pieces occurs on a much shorter 
timescale. 

Paper I suggests the number of orbits required for the as¬ 
teroids that were modeled to dissociate completely into its in¬ 
dividual constituents was tens to hundreds. The pericentres of 
these asteroid orbits were within about 10% of the WD’s Roche 
radius, and the modeled asteroids were strengthless rubble piles. 
If the asteroids were instead internally more cohesive, then the 
disruption limits would change. Asteroids closer to the outer 
edge of the Roche radius may indeed disrupt into sufficiently 
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small bits on longer timescales, such that the Yarkovsky effect 
does play a significant role. 

We discuss this possibility further in Section 6. However, for 
the body of this paper, we continue the story begun in Paper 1 
and consider the subsequent evolution of the rings formed from 
disruption due to close pericentre passages. We assume that the 
constituents of these rings are small enough to be unaffected by 
the Yarkovsky effect, and we operate on this assumption for the 
remainder of Sections 2-5. Consequently, equation © simplifies 
to 

(dv\ Y=0 ALQ , 

\dt ) rad 4-7T mcr 2 

where 

Q = Qabs + Qiet- (4) 

The luminosity of the WD is a steep function of the star’s 
cooling age within the first hundreds of Myr. We adopt the same 
depen dence on cooling age as in equation (6) of lBonsor fe Wvattl 
(l2010lf and re-express their equation as 

£ <‘> = (“- 1 + 5^)"’'“ < 5 ) 

which is applicable for about the first 9 Gyr of WD cooling 
(sufficient for our purposes). Here M* refers to the stellar mass. 

We now determine how the radiative force from equation m 
reshapes the debris rings. First, we consider the perturbations 
on orbital timescales, and then model the cumulative (i.e. long¬ 
term, averaged, or secular) consequences. 


3 UNAVERAGED EQUATIONS OF MOTION 

The orbit of the particle and star is fixed unless acted upon by 
perturbati ons, such a s that from equati o n (l3l) . Based on equa- 
tions from lEfroimskvl d2005l f and iGurfill d2007h , IVeras fe Evansl 
ll2013al l outlined an algorithm which can produce expressions in 
orbital elements for the time evolution of a (semimajor axis), e 
(eccentricity), i (inclination), Q (longitude of ascending node), 
to, (argument of pericentre), / (true anomaly), and the pericen¬ 
tre q = a(l — e) and apocentre q = a(l + e), in the perturbed 
two-body problem. This technique eliminates all Cartesian com¬ 
ponents, and does not make any assumptions about the orbital 
elements except that the motion remains bounded. This facet 
is important because the eccentricities we will be treating are 
near unity, and hence traditional disturbing functions expanded 
about small eccentricities would not be applicable^- 

The complete equations of motion due to th e perturba¬ 
tion fro m equation © have already been derived in IVeras et alj 
(l2015a| j. and we do not repeat them here. One quantity of par¬ 
ticular interest not emphasized in that paper is the closeness of 
the orbital pericentre to the WD for an extremely eccentric or¬ 
bit. We obtain the time evolution of q here by manipul a ting th e 
equations for the evolution of a and e in IVeras et al.l (I2015af l , 
and find 


2 IVerad d2007tf showed how increas ing the order of the expansion helps 
only until the Sundman criterion ilFciTaz-Mclldll 1 994]) is violated. 



ALQ (1 + e cos /) 2 
2-7T m (1 + e) 2 


sin / / 1\ 

2 no? \J\ — e 2 \ c / 


sin 2 (■£) [2 — e — ecos/] 
a (1 — e 2 ) 



( 6 ) 


We need to know how much radiation will drag a parti¬ 
cle located at the pericentre of the orbit (/ = 0) towards (or 
away from) the WD relative to the particle’s unperturbed oscu¬ 
lating location. In the next section we can obtain an order-of- 
magnitude estimate of the maximum value of this displacement 
by averaging equation © over an entire orbit. 


4 AVERAGED EQUATIONS OF MOTION 


In order to consider the cumulative effect of WD radiation over 
many orbits, we analyze the averaged equations of motion. We 
perform the averaging over the true anomaly such that for an 
arbitrary variable /3, 


(dp\ _ J_ [ 2n dp (l-e 2 ) 3/2 
\ dt / 2tv J 0 dt (l + ecos/) 2 

We obtain (IVeras et al.ll2015ah 

/ ( d -iT°\ - ALQ ( 2 +3g2 ) 

\\dtj I 4nmac 2 (1 — e 2 ) 3 ^ 2 


(7) 

( 8 ) 


/(deV= 0 \_ 5ALQe 
\\dt J j 8nma 2 c 2 y/l — e 2 

Equations © and © agree with lWvatt fe Whippld (Il950l h and 
demonstrate that radiation decreases both the semimajor axis 
and eccentricity over time. Note that the (1/c) terms vanish 
upon averaging. 

There is no long-term change in the argument or longitude 
of pericentre 



( 10 ) 


which importantly illustrates that the long-term precession of 
the orbit is dictated by general relativity (in the absence of 
other, additional forces). Finally, we find that 


/ Idq\^ °\ _ ALQyJl — e 2 (4 — e) 
\ V dt ) j 8-irmac 2 (1 + e) 2 

/ (_ ALQV 1 - e 2 (4 + e) 
\\dtj J 8nmac 2 (1 — e) 2 


proving that Poynting-Robertson drag decreases both the peri¬ 
centre and apocentre over time. 

A closer look at equation CD can reveal by how much radi¬ 
ation causes the pericentre to shrink per orbit. Let A represent 
the actual close encounter distance minus the distance predicted 
by Newtonian gravity with no additional forces per orbit. We 
find 

Vl — e 2 (4 — e) \ / at 1 / 2 / m \ 5 

(1 + e) 2 / v 10 au/ \ 10 12 kg / 


A -0.364 m Q 
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( M * \ ^ ( p 

\0.6Mq ) V 2 g/cm 3 


2 

3 


0.1 + 


Myr/ 


(13) 


where p represents the particle density. Hence, along a single 
orbit, radiation from a newly-born WD drags a body towards 
the star by about 4 m at the pericentre. For cooling ages of 1 
Myr and 1 Gyr, this value drops to about 0.4 m and 0.4 mm. For 
highly eccentric orbits, these estimates would be reduced by at 
least one order of magnitude. To place these values in context, 
consider that a typical WD radius is ~ 0.015 - Rq = 10 4 km 
(lHamada fe Salpeterj Il96ll : iHolberg et al.l 120121 ; IParsons et all 
1201211 and the maximum possible extent of the WD disruption 
sphere is ~ 3.6Rq = 2.5 x 10 6 km (Paper I). Hence, over the 
lifetime of the Universe, this gradual pericentre drift is negligible 
compared to the WD radius or its disruption sphere. 

Further, we know that general relativity produces no such 
accumulation of inward or outward drag despite causing a po¬ 
tentially relatively large deviation from the Newtonian orbit at 
each pericentre passage. The maximum value of this deviation 
is about 9 km and the bo dy is pushed towards the 

star only when e < 0.359 (iVerasll2014ll . 

By using both the averaged equations of motion (equations 
m ei and m and a prescription for the luminosity evolution 
(equation 0) , we can numerically determine the long-term or¬ 
bital evolution of particles from WD radiation. This approach 
requires solving coupled differential equations, as we cannot find 
an explicit solution of the general equations. However, in the 
specific case of high eccentricity, we have found an explicit so¬ 
lution, as described in the next section. 

First, we numerically integrate these full averaged equations 
of motion (equations O El El and 1101) to determine the evolution 
for WD planetary systems. Results of this integration are plot¬ 
ted as solid lines in Figure [I] for bodies with radii ( R ) spanning 
four orders of magnitude. Only for this size range (10 -5 — 10 _1 
m) can we be sure that other radiative effects do not come into 
play. Each panel illustrates a different WD cooling age ti n i; older 
WDs take longer to shrink orbits. For simplicity, we adopted 
a characteristic WD mass (M*) of 0.6 Mq, characteristic body 
density (p) of 2 g/cm 3 and Q = 1, for the integrations. As¬ 
suming Q = 1 means that no directed scattering or reemission 
occurs (i.e. Q le f = Q yar = 0), such that all radiation is absorbed 
and emitted omnidirectionally. This assumption accurately re¬ 
produces the properties of small carbonaceous dust particles 
with albedos of about 0.005. Alternatively, typical asteroids with 
albedos between 0.1 and 0.3 would satisfy 1.1 < Q < 1.30 

We also adopted one of the pericentre values used in Paper 
I to better demonstrate the link with that paper and to pro¬ 
vide analytically-motivated initial parameters. In that respect, 
q m in = 0.0126-Rq = 5.86 x 10~ 5 au is the radius of a typical 
0.6 Mq WD, and g max = 2.73 Rq = 0.013 au is the maximum 
value of the disruption sphere for that WD mass. 

The figure demonstrates that for a wide range of WD cool¬ 
ing ages, all particles in this size range will maintain their orig¬ 
inal semimajor axes for the vast majority of their evolution, 
before their orbits suddenly shrink to within the WD Roche ra¬ 
dius. Also plotted is the apocentre of the orbit. For the majority 


3 The albedo distribution for main belt asteroids is bimodal, with 
peaks at about 0.05 and 0.25. 


of the particle’s evolution, the particle maintains its initial ec¬ 
centricity; circularization occurs only as the semimajor axis is 
shrinking drastically. This behaviour is mirrored for particle evo¬ 
lution at later WD cooling ages, although the collision timescale 
is increased, as will be shown in the next section. 


5 HIGH ECCENTRICITY LIMIT 


Numerical simulations from the last subsection demonstrate 
that a particle on an initially highly eccentric orbit will remain 
so for the vast majority of its lifetime. This result suggests that 
considering the high-eccentricity limit of the equations might 
represent a good approximation to the actual motion. Our ob¬ 
jective in this subsection is to express this motion by an explicit 
algebraic expression, which will facilitate future study. 

We perform series expansions about t = 1 — e, where e< 1, 
and denote the result with a tilde. We obtain 


5 ALQ 


8y/2nmac 2 (1 — e) 3 ^ 2 


9ALQ 


32\/27r mac 2 (1 — e) 1 ^ 2 


+ o((l-e) 1/2 ), (14) 


deV 
dt) 


5 ALQ 


8V2nma 2 c 2 (1 — e) 1 ^ 2 


+ C>((l-e) 1/2 ). 


(15) 


Now suppose the debris ring is formed at the WD cooling age 
tini, such that at this time each particle adopts semimajor axis 
and eccentricity values of ai n i and e ln i. We see from equation 
(1111) that 



(16) 


The semimajor axis and eccentricity change due to radiative 
effects, but on average the pericentre remains constant. The 
eccentricity and semimajor axis evolution can thus be approxi¬ 
mated as 


aft) [1 - e(t)] ^ a in i (1 - e in i) • (17) 

Equation (El may also be obtained by simultaneously solving 
the leading terms present in equations El and EI- By using 
the relation in equation El- we obtain 


/ (daV °\ __ 5 ALQVa _ 

\\dt J J 8 V^nmaf^ c 2 (1 - e in i) 3/2 

In order to solve equation El algebraically and explicitly 
for the semimajor axis, we express L as in equation © but 
replace the power law exponent of —1.18 with —6/5, such that 


L 


Lo 


± M" 6/5 

10 + Myr J 


(19) 


where Lo = 3.26Lq(M/0.6Mq) and is hence constant. Also, we 
express A and m in terms of the body’s radius (R) and density 
(p), as we are more interested in the dependence on size than 
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Orbit shrinkage from radiation 

Qini = Qmin't‘'* /o (q max ~ Q m i n) 




Figure 1. Orbit shrinkage of various particle sizes (with radii R) 
from WD radiation. The three panels refer to different WD ages 
(iini) for a WD with a typical mass of 0.6Mq. From top to bottom, 
these ages correspond to luminosities of about 22 Lq, 2 x 10“ 3 Lq and 
1 x 1O _4 L 0 . Plotted in each pair of similarly-coloured curves are the 
evolution of the semimajor axis (a; lower curve) and apocentre ( q ; 
upper curve). The eccentricity of the orbit becomes zero only when 
those two lines converge. Solid lines are the result of numerical in¬ 
tegrations of the full equation of motion m and the dashed lines 
are from our analytical approximations (equations I21I22D . The ini¬ 
tial semimajor axis was taken to be ai n i = 10 au and the pericentre 
to be < 7 i n i = q m [ n + 5%(q , m ax — 9min), which corresponds to the or¬ 
bital parameters of some rings formed in Paper I. All orbits become 
compressed to entirely within the Roche radius. 


on mass. The final result can be expressed compactly through 
the following auxiliary variable 


f 75-5 1/s \ f LoQMyr 

\ 32 ■ 2 3 / 10 7t/ {<?pRa ini ( 1 - e ini ) 3 / 2 


X 



10tuu\ 1/5 

Myr ) 




The evolution of the semimajor axis is 


( 20 ) 


xr 2 

a(t) = Oini — K + - -. (21) 

dUini 

Further, because the pericentre is constant, the apocentre reads 


q(t) a(t) - 2<?ini 


( 22 ) 


We plot the curves predicted from equation (EH as dashed 
lines on Fig.Q] The agreement with the true solution is excellent, 
for both the semimajor axis and apocentre, and for all cooling 
ages relevant to WD pollution, at least to the precision we seek. 

Only physical solutions from equation EH) are plotted. Af¬ 
ter a body encounters the Roche radius, the time evolution is 
stopped. Otherwise, the semimajor axis would then unphysically 
increase after achieving a minimum. To stop the time evolution 
at the appropriate time, one can multiply equation EH by the 
appropriate Heaviside function. 

Of particular interest is the shrinking time , t s h r , which we 
define as the time taken for the particle’s initial orbit to be com¬ 
pressed entirely within the WD Roche radius. By using equation 
(1211) . we find 

/ 10(t s hr+iini)V V5 A 10tini\" 1/5 
V Myr ) V Myr ) 


/ 64 • 2 3/10 tt\ 
V 75-5 1 / 5 ) 



ini + 7?Roche) (23) 


where RR OC he is the Roche radius of the WD. From this equa¬ 
tion fshr may be solved for explicitly; we adopted the (physical) 
solution with the lower sign. 

We use equation (1231) to generate Figure [H which illustrates 
the dependencies of t s hr on R and a. For observational perspec¬ 
tive we display the range of coo ling ages of known WD with de- 
bris discs (as 30 Myr - 1.5 Gvr: iFarihi et al.ll2009l ; iGirven et ahl 
l2f)12l : lBergfors et al.ll2fll3 l with a thick brown horizontal bar and 
vertical brown lines. The steep increase in shrinking timescale 
during the first Gyr of WD evolution is a reflection of the star’s 
rapidly dwindling luminosity at early cooling ages. Neverthe¬ 
less, for all cooling ages, the radiation is still strong enough to 
shrink and circularize the orbits on timescales orders of mag¬ 
nitude smaller than the WD cooling age. Coincidentally, the 
largest particles for which the Yarkovsky effect play no role 
(~ 10“ 1 m) also approach the maximum particle size at which 
f s hr <S fini, at least for ai n ; = 10 au. Indeed, we do not yet have 
any observational constraints on ai n ;; equations such as EH ma y 
help motivate likely values. 


6 DISCUSSION 

With our results, we can now try to construct a consistent pic¬ 
ture from Paper I about how tidal disruption of asteroids is 
followed by orbit compression. A single asteroid is tidally dis¬ 
rupted into a collisionless ring with an unknown size distribu¬ 
tion. The orbit of the ring is equivalent to the orbit of the pro¬ 
genitor, whose eccentricity exceeds 0.98 and whose semimajor 
axis exceeds a few au. A Main Belt-like asteroid would satisfy 
aini ~ 5 au, a Kuiper Belt-like asteroid would satisfy a ini « 30 
au, and an asteroid from a debris belt in a system like Fomal- 
haut would satisfy a; n ; ~ 100 au. Because the ring is assumed to 
be collisionless, the bodies will not gravitationally perturb one 
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1<T 3 1(T 2 10- 1 10° 10 1 10 2 10 3 10 4 

tini / My r 


Figure 2. The shrinking timescale t s h r as a function of WD cooling 
age ti n i f° r curves of different particle radii (upper panel) and initial 
semimajor axes (lower panel). In all cases, a p = 2 g/cm 3 particle 
orbits a 0.6Mq WD. In the upper panel, we adopt ai n i = 10 au 
and < 7 i n i = g m i n + 10%(<? ma x - 9min) (solid lines) and <? in i = q min + 
5%(^max — 9min) (dashed lines). The similarly-coloured pairs of lines 
each corresponds to a different order of magnitude of particle radius. 
In the lower panel, all particles have R = 10~ 3 m and q; n ; = q m ; n + 
5%(Y/n, flx — (Jmin)- Marked for reference as a thick brown bar and 
vertical brown lines is the range of cooling ages at which WD discs 
have been observed. All curves in both panels demonstrate that the 
shrinking timescales are orders of magnitude shorter than the cooling 
times. 


another. Every body will experience radiation forces in the form 
of Poynting-Robertson drag. 


6.1 Collisions upon contraction 

If the ring contains differently-sized bodies, then the effects of 
Poynting-Robertson-drag alone will cause the orbits of these 
bodies to contract at different rates. As smaller bodies contract 
more quickly, their pericentres will precess more rapidly due to 
general relativity. Consequently, collisions might occur. Consider 
equation 19 of Paper I: 


Pu, 


0.15 Myr 


1 - e 2 
1 - 0.999 2 


f Mwd \ ~ 3/2 / a \ 5/2 
\0.6 Mq J \ 1 au/ 


(24) 


which illustrates that the general relativistic precession period 
Pu, oc a 5 / 2 and is also a function of e and Mwd- Both these last 
two quantities remain nearly constant as the orbit contracts by 
a factor of tens to hundred^} Therefore, suppose an asteroid on 
a a = 10 au orbit is broken up into large and small particles with 
radii Ri alge and i? sma n. The larger particle ring will continue to 
precess with a period of about 47 Myr while the smaller particle 
ring is shrunk; at a = 1 au, that small particle ring’s precession 
period is just about 0.15 Myr. Consequently, when t s hr > Pu,, we 
can expect collisions to be significant; for small enough particles 
such that t s hr -C Pu ,, collisional effects would be negligible. 


6.2 Yarkovsky perturbations on chunky rings 


If there exist ring constituents that remain greater than about 
one metre in diameter over long-enough timescales, then the 
Yarkovsky effect will likely become the dominant contribution 
to the motion. We mentioned this possibility in Section 2, but 
now provide some quantification. 

Let Q = Qabsl + Qrefl + kQ yal Y . The components of Q are 
a function of spin state and orbital position, velocity and time. 
Al though the s implified spherical internal heat model provided 
bv iBro 1 (l2006tl allows the time evolution of the semimajor axis 
or eccentr icit y to be expressed ( see Sections 2.3.2 and 2.3.3 of 
lBrozll2006l. or Vo krouhlickvll998l . lVokrouhlickv fc Farinellal[l998l 


or Vokrouh lickv 199^), that model itself might not be applicable 

for the extremely high eccentricities considered in this paper. For 
a particle on an eccentric orbit, the seasonal component of the 
Yarkovsky effect contains many forcing periods; the amplitudes 
of those frequencies are functions of the orbital eccentricity. 

Regardless of these difficulties, we now provide a zeroth- 
order demonstration of the potential importance of the effect. 
Each element of Q is strictly bounded between 0 and 2, allow¬ 
ing us to bound the extent of the potential perturbation. By 
assuming th at the matrix elemen ts are independent of position 
and velocity, IVeras et al.l (l2015al l obtained the resulting orbital 
element-based unaveraged equations, as well as the averaged 
equations to leading order in (1/c). By expanding their equa¬ 
tions A3 and A4 about e = 1 - e, £ C 1, we obtain 


AL 


Ancmna 2 


1 | 1 

2(1 —e) + 4 


cos* 

sin * sin Q, 
sin * cos O 


+ 0(1 —e), 


(25) 


AL 


8ncmna 3 


cos* 

sin * sin Q, 
sin * cos Q 


+ O (Vl^e) ■ (26) 


Comparing equations (1251) - (12611 with (I14l) - (115l) reveals that (1141) - 
(1151) are a factor of about 1/c smaller than (125ll - (l26ll . and can 


4 Eventually, the eccentricity starts to decrease appreciably, poten- 
tially increasing the factor in square brackets by up to a factor of 
about 500. Nevertheless, this eccentricity decrease does not occur un¬ 
til the particle has migrated nearly all the way to the Roche radius. 
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be neglected when the Yarkovsky acceleration is “turned on”. 
Further, because the leading-order eccentricity term of equation 
1(25)1 is of order 1/(1 — e) whereas that of equation (1261) is of 
order of just unity, the Yarkovsky acceleration might cause the 
orbit to shrink quickly at a nearly fixed eccentricity. 

The off-diagonal terms will determine the sign of the expres¬ 
sions, and hence whether the orbit shrinks or expands, circular¬ 
izes or becomes parabolic. In fact, because the initial eccentricity 
is so high, on average we might expect about half of the large 
chunks in the ring to be ejected. Given that the Yarkovsky ac¬ 
celeration affects particles only above a certain size, the extra 
movement of these particles with respect to the relative velocity 
of the constituents after break-up might lead to collisions, which 
could produce a damping effect. 

Although we wish to obtain an expression similar to equa¬ 
tion for the Yarkovsky effect, we cannot do so, but can 
make some progress towards a solution. In principle, we could 
follow the same procedure as in Section 5 to obtain an approx¬ 
imate closed formula for a{t) due to the leading-order terms in 
equations (1251) - 11261) . In the high eccentricity limit, we can relate 
a and e to one another and their values at a particular cooling 
age because of the similar forms of equations (1251) - (1261) . Dividing 
equation (HD by equation (EH) and integrating yields 

2a 2 = 2a 2 ni + ± (e - e ini ) + In ) ' (27) 

However, we cannot go further and substitute equation (1271) into 
equation (1251) and then integrate primarily because equations 
dlOi-Ei) depend on the inclination and longitude of ascending 
node (which are time-dependent) and the elements of Q (which 
may be time-dependent). 


7 CONCLUSIONS 

We have addressed a problem in the formation of discs which 
orbit WDs at a distance of IRq: how to shrink the 10 — 100 
au-scale extremely eccentric (e > 0.98) bound orbits of debris 
from tidally-disrupted asteroids. We demonstrated that this or¬ 
bit compression readily occurs for particle sizes of 10~ 5 — 10 _1 
m on timescales many orders of magnitude shorter than the WD 
cooling age, due to WD radiation alone. We provided explicit ap¬ 
proximations for the shrinking time ('equation 1231) as well as for 
the semimajor axis and apocentre evolution (equations 12111221) 
during the orbit compression (the pericentre remains constant). 
After the debris is perturbed close to or within the Roche ra¬ 
dius, subsequent evolution is likely to be dictated by collisional 
and sublimative forces. 
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